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Flavor oscillations of neutral B mesons have been studied in e'^'e” annihilation data collected 
with the BABAR detector at center-of-mass energies near the T(4S') resonance. The data sample 
used for this purpose consists of events in which one B" meson is reconstructed in a hadronic decay 
mode, while the flavor of the recoiling B" is determined with a tagging algorithm that exploits 
the correlation between the flavor of the heavy quark and the charges of its decay products. From 
the time development of the observed mixed and unmixed final states we determine the B^-B” 
oscillation frequency Am^ to be 0.516 ± 0.016(stat) ± O.OlO(syst) ps“^. 

PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 

In the Standard Model, B^-B” mixing Q occurs change of up-type quarks, with the top quark contribut- 

through second-order weak diagrams involving the ex- ing the dominant amplitude. A measurement of the mass 
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difference Am^ between the mass eigenstates is there¬ 
fore sensitive to the value of the Cabibbo-Kobayashi- 
Maskawa [Q matrix element Vtd- Tbe phenomenon of 
particle-anti-particle mixing in the neutral B meson sys¬ 
tem was first seen almost fifteen years ago |^. The os¬ 
cillation frequency Am^ has been measured with both 
time-integrated and time-dependent techniques Q. 

In this Letter we present a measurement of time- 
dependent mixing based on a sample of 32 million BB 
pairs recorded at the T’(4iS') resonance with the BABAR 
detector at the Stanford Linear Accelerator Center. This 
study and a related CP asymmetry measurement [|j are 
described in more detail in Ref. [^. At the PEP-II 
asymmetric-energy e+e“ collider, the T(4S') provides a 
source of B^B^ pairs moving along the e~ beam direc¬ 
tion (z axis) with a known Lorentz boost of /Jy = 0.55, 
allowing a novel technique for determining Anid- 

The B^-B^ mixing probability, for a given B^ lifetime 
r, is a function of Am^ and the proper decay-time differ¬ 
ence At between the two neutral B mesons produced in 
a coherent P-wave state in the T(4S') event. The result 
is a time-dependent probability to observe unmixed (-I-), 
B^B^, or mixed (—), B^B^ and P°P°, events: 

Prob{B°B° B°B°,B°B° or WB°) oc 

e-l^*l/’^(l ± cos AmdAt). (1) 

The effect can be measured by reconstructing one B 
in a flavor eigenstate, referred to as Prec, while the 
remaining charged particles originating from the decay 
of the other B, referred to as Ptag, are used to iden¬ 
tify, or “tag”, its flavor as a P° or B^. The charges 
of identified leptons and kaons are the primary indica¬ 
tors, although other information in the event can also be 
used to identify the flavor of Ptag- The time difference 
At = Pec —ttag — Az/Pjcis determined from the separa¬ 
tion Az of the decay vertices for the flavor-eigenstate and 
tagging P decays along the boost direction. The tagging 
and vertexing algorithms used in this analysis are nearly 
identical to those employed for CP violation studies, in 
which one P is fully reconstructed in a CP eigenstate [|j . 

The value of Amu is extracted from a tagged flavor- 
eigenstate B^ sample with a simultaneous unbinned max¬ 
imum likelihood fit to the At distributions of mixed and 
unmixed events. There are two principal experimental 
factors that complicate the probability distribution given 
by Eq. (|^. First, the tagging algorithm, which classifies 
events into categories i depending on the source of the 
available tagging information, incorrectly identifies the 
flavor of Ptag with a probability Wi. This mistag rate 
reduces the observed amplitude of the oscillation by a 
factor (1 — 2wi). Second, the resolution for At is com¬ 
parable to the oscillation period and must be well un¬ 
derstood. The probability density functions (PDFs) for 
the unmixed and mixed signal events, P±,sig, can be ex¬ 
pressed as the convolution of the underlying At distribu¬ 


tion for the tagging category, 

g-|At|/r 

h±{At; Amd, Wi) = —--[1 ± (1 — 2wi) cos ArudAt] , 

4r 

with a At resolution function containing parameters dj. 
A log-likelihood function is then constructed by summing 
lnP±,sig over all events within each of the tagging cat¬ 
egories. The likelihood is maximized to extract simul¬ 
taneously the mistag rates Wi, the resolution function 
parameters dj, and the mixing parameter Arud- 

The BABAR detector is described in detail elsewhere . 
Charged particles are detected and their momenta mea¬ 
sured by a combination of a 40-layer drift chamber 
(DCH) and a five-layer silicon vertex tracker (SVT) em¬ 
bedded in a 1.5-T solenoidal magnetic field. Prec de¬ 
cay vertices are reconstructed with a resolution of typ¬ 
ically 65 fj,m along the boost direction. A detector of 
internally reflected Cherenkov radiation (DIRC) is used 
for charged hadron identification. Kaons are identified 
with a neural network based on the likelihood ratios in 
the SVT and DCH, derived from dE/dx measurements, 
and in the DIRC, calculated by comparing the observed 
and expected pattern of Cherenkov light for either kaons 
or pions. A finely segmented CsI(Tl) electromagnetic 
calorimeter (EMC) is used to detect photons and neu¬ 
tral hadrons, and to identify electrons. Electron can¬ 
didates are required to have a ratio of EMC energy to 
track momentum, an EMC cluster shape, DCH dE/dx, 
and DIRC Cherenkov angle consistent with expectation. 
The instrumented flux return (IFR) is segmented and 
contains resistive plate chambers for muon and neutral 
hadron identification. Muon candidates are required to 
have IFR hits located along the extrapolated DCH track, 
an IFR penetration length, and an energy deposit in the 
EMC consistent with the muon hypothesis. 

Neutral P mesons are reconstructed in a sample of 
multihadron events in the flavor eigenstate decay modes 
£)(*)- 7 r+, and J/ipK*^. The decay 

channels AT+tt”, AT+tt^tt^, Ar+7r+7r“7r“ and AlgTr+Tr” 
are used to reconstruct P° candidates, while the modes 
Ar+7r“7r“ and AT^tt" are used for D~ candidates. 
Charged D*~ candidates are formed by combining a 
P° with a soft 7r“. Finally, the P° candidates are 
formed by combining a D*~ or D~ candidate with a 
7r+, (p+ ^ 7r+7r°) or a)'’ (a)^ ^ 7r+7r“7r+); likewise, 

P° ^ J/ip K*^ candidates are reconstructed from combi¬ 
nations of J/ip candidates, in the decay modes e“'"e“ and 
pC, with a K*^ (Ar*° ^ AT+tt”). The selection and 
reconstruction of these decay chains and the selection of 
multihadron events, including continuum supression, is 
described in more detail in Ref. [|j. 

Neutral P candidates are identified by the difference 
AE between the energy of the candidate and the beam 
energy \/s/2 in the center-of-mass frame, and the beam- 
energy substituted mass tties, calculated from -fyi/2 and 
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the reconstructed momentum of the B candidate. Can¬ 
didates are selected by requiring toes > 5.2 GeV/c^ 
and AE within ±2.5 standard deviations of 0 (typically 
|Ai?| < 40MeV). When multiple candidates in a given 
event are selected (with probability of about 0.25%), only 
the one with the smallest |Ai?| is retained. 

After the daughter tracks of the Brec are removed, the 
remaining tracks are analyzed to determine the flavor 
of the Btag. For this purpose, we use the flavor tagging 
information carried by primary leptons from semileptonic 
B decays, charged kaons, soft pions from D* decays, and 
high momentum charged particles, to assign an event to 
a single tagging category. 

Events are assigned a Lepton tag if they contain 
an identified lepton with a center-of-mass momentum 
greater than 1.0 or l.lGeV/c for electrons and muons, 
respectively. The momentum requirement selects mostly 
primary leptons by suppressing opposite-sign leptons 
from semileptonic charm decays. If the sum of charges 
of all identified kaons is non-zero, the event is assigned 
a Kaon tag. The final two tags involve a multivariable 
analysis based on a neural network, which is trained to 
identify primary leptons, kaons, and soft pions, and the 
momentum and charge of the track with the maximum 
center-of-mass momentum. Depending on the output of 
the neural net, events are assigned either an NTl (more 
certain) tag, an NT2 (less certain) tag, or are not tagged 
at all (about 30% of all events) and excluded from the 
analysis. 

Tagging assignments are made mutually exclusive by 
the hierarchical use of the tags. Events with a Lepton tag 
and no conflicting Kaon tag are assigned to the Lepton 
category. If no Lepton tag exists, but the event has a 
Kaon tag, it is assigned to the Kaon category. Otherwise 
events with neural network tags are assigned to corre¬ 
sponding neural network categories. 

The decay time difference At between B decays is de¬ 
termined from the measured separation Az = tree ~ -^tag 
along the z axis between the reconstructed Brec(^rec) and 
flavor-tagging decay Btag (.Ztag) vertex. This measured 
Az is converted into At with the use of the known 7’(4S') 
boost, including a correction on an event-by-event basis 
for the direction of the B mesons with respect to the 
z direction in the T’(4S') frame. The At resolution is 
dominated by the z resolution of the tag vertex position. 
Reconstruction of the Btag decay vertex starts with all 
tracks in the event except those incorporated in Brec- In 
order to reduce the contamination from D meson decay 
products, those identified as kaons are also excluded. An 
additional constraint is provided by the calculated Btag 
production point and three-momentum, determined from 
the momentum of the B^ec candidate, its decay vertex, 
the average position of the interaction point, and the 
T(4S') boost. Tracks with a large contribution to the 
are iteratively removed from the fit until those remain¬ 
ing (> 1) have a reasonable fit probability or all tracks 


are removed. Only events with a reconstructed Btag ver¬ 
tex, |At| < 20 ps and tXAt < 1.4 ps are retained (about 
84%), where dAt is the measurement error derived from 
the vertex fits. 



FIG. 1: Distribution of mss for all B° candidates with a flavor 
tag and a reconstructed tag vertex. 

The distribution of toes for the selected candidates is 
shown in Fig. |^, where the result of a fit with a Gaussian 
distribution for the signal and an ARGUS function |^] 
for the background is also displayed. The fitted number 
of signal events and their purity (for toes > 5.27 GeV/c^) 
are 6347±89 and (85.8±0.5)%, respectively. The sample 
composition by tagging category is given in Table 


TABLE I: Signal yields per tagging category, obtained from 
the Toes distributions after all selection requirements. The 
purity is quoted for toes > 5.27 GeV/c^. 


Category 

Tagged 

Purity (%) 

Lepton 

1097 ± 34 

96.0 ±0.7 

Kaon 

3156 ±63 

84.6 ±0.7 

NTl 

798 ± 31 

88.9 ± 1.2 

NT2 

1293 ±43 

79.4 ± 1.3 


In the likelihood fit, the At resolution function is ap¬ 
proximated by a sum of three Gaussian distributions 
(core, tail, and outlier) with different means and widths, 


77(<5t;a) = ^ 


fk 


‘5'fcrrAi'\/27r 

fs 


(T3v/27r 


exp 


exp - 




((5t - bkaAtf 
2{Sk<JAt)^ 


2(73^ 


where St — At — Attrue- The sum of the fractions fk is 
constrained to unity. For the core and tail Gaussians, the 
widths Ok = Sk y. <JAt are the event-by-event measure¬ 
ment errors scaled by a common factor Sk- The scale 
factor of the tail Gaussian is fixed to the Monte Garlo 
value since it is strongly correlated with the other res¬ 
olution function parameters. The third Gaussian, with 
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a fixed width of 0-3 = 8 ps, accounts for outlier events 
with incorrectly reconstructed vertices (less than 1 % of 
events). A separate core bias coefficient bi i is allowed for 
each tagging category i to account for small shifts due to 
inclusion of charm decay products in the tag vertex, while 
a common bias coefficient &2 is used for the tail compo¬ 
nent. These offsets are proportional to a At since both 
the size of the bias and the resolution for ztag depend 
kinematically on the polar angle of the flight direction of 
the charm daughter. The tail and outlier fractions and 
the scale factors are assumed to be the same for all decay 
modes, since the Ztag measurement dominates the reso¬ 
lution for At. This assumption is confirmed by Monte 
Carlo studies. Separate resolution parameters are used 
for two different data-reconstruction periods, referred to 
as Runl and Run2, which mainly differ in vertex perfor¬ 
mance and tracking efficiency. 

In the presence of backgrounds, additional terms must 
be added to the signal PDF H±,sig for each background 
source, 

— fi,sigT~(-±,sig T ^ j), 

j=bkgd 

where the background PDFs B±^ij provide an empiri¬ 
cal description for the At distribution of the background 
events in each tagging category i of the sample. The frac¬ 
tion of background events for each tagging category and 
background source is given by /ij, while are pa¬ 

rameters used to characterize each source of background 
by tagging category for mixed and unmixed events. The 
signal probability /i,sig is determined from the measured 
event mss on the basis of a separate fit to the observed 
Toes distribution in tagging category i. The sum of signal 
and background fractions is forced to unity. 

The At distributions of the background are described 
with a zero lifetime component and a non-oscillatory 
component with non-zero lifetime. We fit for separate 
resolution function parameters for the signal and the 
background to minimize correlations. The likelihood fit 
involves a total of 44 parameters, including Am^, the 
average mistag fraction and the difference between 
and mistags for each tagging category ( 8 ), parame¬ 
ters for the signal At resolution (16), and parameters for 
background time dependence (5), At resolution ( 6 ), and 
effective dilutions ( 8 ). The value of Am^ was kept hid¬ 
den throughout the analysis until all analysis details and 
the systematic errors were finalized, to eliminate possible 
experimenter’s bias. 

The results from the likelihood fit to the tagged B^ 
sample are summarized in Table |I[ The probability to 
obtain a likelihood smaller than that observed is 44%, 
evaluated with a parameterized Monte Carlo technique. 
The value of Am^ given by the fit, prior to final correc¬ 
tions, is Arridfit = 0.525 ± 0.016 ps“^. One method for 
displaying the result of the full likelihood fit is to use the 


TABLE II: Results for Arud and a subset of the parameters 
obtained from the likelihood fit to the At distributions. ArUd 
includes small corrections described in the text. 


Parameter 

Fit Value 

Parameter Fit Value 

Anid 

0.516 

± 

0.016 





w (Lepton) 

0.079 

± 

0.014 

w (NTl) 

0.219 

± 

0.022 

w (Kaon) 

0.166 

± 

0.012 

w (NT2) 

0.344 

± 

0.020 

5i(Runl) 

1.37 

± 

0.09 

5i(Run2) 

1.19 

± 

0.11 

/ 2 (Runl) 

0.014 

± 

0.010 

/2 (Run2) 

0.015 

± 

0.010 

/a (Runl) 

0.008 

± 

0.004 

/3(Run2) 

0.000 

± 

0.014 




lAtl (ps) 

FIG. 2: Distributions of At in data for the selected a) un¬ 
mixed and b) mixed events (mEs(Brec) > 5.27GeV/c^), with 
projections of the likelihood fit (solid) and the contribution of 
the background (dashed) overlaid. The time-dependent mix¬ 
ing asymmetry Amix{\At\) is shown in c). 


observed mixing asymmetry, 

^ (At) JV^txedi^At') 

— — TXB' 

unmixedy^^) “r lyrnixedy^^) 

If the flavor tagging and At determination were perfect, 
the asymmetry as a function of At would be a cosine with 
unit amplitude. The amplitude is diluted by the mistag 
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probability and the experimental resolution for At. The 
observed At distributions of both the mixed and unmixed 
events, and their asymmetry Amix are shown along with 
projections of the likelihood fit result in Fig. ^ 

Since the parameters of the At resolution for both sig¬ 
nal and backgrounds are free parameters in the fit, their 
contribution to the uncertainty on Arrid is included as 
part of the statistical error. Remaining systematic er¬ 
rors arise from the choice of the signal At resolution de¬ 
scription, its capability to handle outliers and various 
worst-case SVT misalignment scenarios (±0.005ps“^), 
and by approximations and uncertainties in the Az to 
At conversion due to the knowledge of the absolute z 
scale of the BABAR detector and PEP-II boost (less than 
±0.002ps“^). Systematic errors due to background in¬ 
clude the choice of its At distribution and resolution de¬ 
scription (±0.002ps“^), variation of the sum of back¬ 
ground fractions from the separate toes fitSj and the un¬ 
certainty on the magnitude of the small B'^ component 
of the signal (±0.002 ps“^). A correction of—0.002 ps“^, 
derived from data, is made to account for the small varia¬ 
tion of the background composition as a function of toes j 
which affects the background At distribution. The sta¬ 
tistical error (±0.002 ps“^) on this extrapolation from 
the sideband to the signal region is included as a system¬ 
atic uncertainty. An additional correction of —0.007 ps“^ 
is applied for a bias observed in fully simulated Monte 
Carlo events. The bias is mainly due to correlations be¬ 
tween the mistag rate and the At resolution that are 
not explicitly incorporated into the likelihood function. 
The systematic error assigned to this correction includes 
contributions from the statistical precision of the Monte 
Carlo study (±0.003ps“^), model variations due to un¬ 
certain branching fractions and lifetimes of the tag-side 
D mesons and the assumed fraction of wrong-sign kaons 
produced in B decays (±0.001 ps“^), and variation of 
the requirement on the maximum allowed value of UAt 
(±0.003 ps“^). Finally, the variation of the fixed B^ life¬ 
time within the known errors [Q leads to a systematic 
uncertainty of ±0.006 ps“^. 

In conclusion, a new technique involving the time- 
difference distribution of a tagged sample of fully- 
reconstructed neutral B decays has been used to deter¬ 
mine the B^-B^ mixing frequency Ato^ to be 

Arud = 0.516 ± 0.016(stot) ± O.OlO(syst) ps“^. 


This is one of the single most precise measurements avail¬ 
able. Moreover, the error on Ato^ is still dominated 
by the size of the reconstructed B^ sample, leaving sub¬ 
stantial room for further improvement as more data are 
accumulated. The result is consistent with the current 
world average 0 and a recent BABAR measurement with 
a dilepton sample |^. The analysis shares the same 
flavor-eigenstate sample as used for the determination of 
sin2/3, thereby providing an essential validation for the 
reported sin2/3 result ||]. 
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